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ABSTRACT 


Real-time detection of R peaks in QRS complex of ECG signal is the first 
step in the processing of ECG waveform. Based on this, various other ECG 
parameters can be extracted. These parameters provide substantial 
information about various heart diseases. In this paper, we are proposing a 
method to detect R-peaks of ECG signal dynamically. The most prominent 
role in the R—peak detector is executed by the microcontroller. This method 
originates by acquiring signal from the subject and necessary pre-processing 
is carried out on the signal in order to achieve the denoised signal. 
Subsequently, this filtered signal is handed over to microcontroller where a 
pulse is generated for each R—peak that is found in the QRS complex of ECG 
signal. The microcontroller is embedded with a signal processing algorithm. 
The algorithm used to determine the R-—peaks is double differentiation 
method which is straightforward and robust. 
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1, INTRODUCTION 

The electrical conductivity of heart is nothing but ECG which can be acquired by placing the 
electrodes in a specific configuration on the surface of the skin. ECG signal gives the vital information about 
the functionality of heart. R wave has its own significance in the QRS complex. By detecting R peaks, the 
first crucial step is finished in analyzing the ECG. The other features are obtained with reference to R-wave. 
The acquired ECG signal on the skin surface is poor in amplitude which 1s around 0.05mv to 3 mV and the 
frequency range of the ECG signal is 0.01 to 300 Hz [1]. In consideration of the weak and noisy ECG signal, 
it has to be boosted without amplifying the noise. Hence, we have used low pass, high pass and notch filters. 
Subsequently, peaks of the ECG signal are to be detected dynamically, using a microcontroller. 

Many methods such as Derivative approach, Wavelet approach, neural network approach, Hilbert 
transform approach are available in the literature [2] for detection of QRS complex. All of these algorithms 
are highly sensitive to noise. Moreover, these algorithms are complex, time-consuming and are difficult to 
implement dynamically on the microcontroller [3-10]. 

In this paper, we have proposed an algorithm based on Double differentiation approach. In this 
process, analogue ECG signal is given to ADC channel of microcontroller where the signal is sampled. This 
signal is nothing nut the discrete signal. Later, first and second differentiation procedures are performed on 
the discrete signal. Next, summing and scaling of the derivatives are carried out for the discrete signals. Next, 
the obtained values are compared with a threshold, 1.e. obtained during the initial phases of processing. This 
threshold value is dynamically derived, as the signal characteristics vary from subject to subject. This process 
provides an output on the I/O pin of a microcontroller as and when the R-wave appears in the ECG signal. 
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2. DESIGN AND METHODOLOGY 
2.1. Design 

In this proposed work, we have designed an ECG amplifier for acquiring and pre-processing the 
signal. The block diagram of ECG amplifier is shown in Figure 1. 

In Figure 1, it is observed that the ECG signal is acquired from the subject’s body by connecting 
electrodes in lead II configuration. The acquired signal from the subject is very low in amplitude which is 
around 0.05 to SmV and easily distorted by movement artefacts and other noise sources like power lines. In 
this procedure, retaining the signal without distortion is important. In consequence, the signal is amplified by 
the instrumentation amplifier whose gain is large enough to amplify the signal. This signal is then passed 
through high pass filter with a cut-off frequency of 0.5 Hz to remove the noise caused by electrode half-cell 
potential variations. 

Subsequently, the signal is passed through notch filter to remove the mains noise whose frequency is 
50Hz. Once again, the signal is passed through low pass filter with a cut-off frequency of 1 KHz to remove 
redundant high frequency components that are present in the signal. Finally, the denoised signal is obtained. 
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Figure 1. ECG Amplifier 


Considering that, we are feeding the ECG signal to ADC channel of the microcontroller. The input 
must be within the ADC limits and it should be unipolar in nature. Here, an amplifier is used to raise the 
signal amplitude so that the signal amplitude is within+2.5 V. Then the signal is fed through a level shifter, 
so that the signal is within 0-5V range. The output of the ECG amplifier can be seen in Figure 2. 
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Figure 2. Signal obtained from ECG Amplifier. 


The microcontroller in the R—peak detector is programmed with signal processing algorithm to 
detect the R—peaks of ECG signal. The ADC channel of the microcontroller samples the analog signal and 
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with necessary signal processing techniques, a pulse is generated at each R—peak in QRS complex of ECG 
signal which is shown in Figure 3. 
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Figure 3. Detection of R-peaks. 


2.2. Methodology 

In Biomedical signal processing text book [11], the author proposed a method on double 
differentiation. We have modified that algorithm in order to perform a dynamic operation on the discrete 
form of ECG signal. The proposed algorithm is simple to implement. This process involves three stages 
which are explained in further sections. 


2.1.1 Shifting and Smoothing of the data: 

For every 2ms, interrupt is generated from the built-in timer to receive the sample from ADC 
channel which is stored in X [OQ]. In the process of shifting of data, we have collected five samples from ADC 
channel. Though we have used filters in the hardware section, because of differentiation, high frequency 
noise gets amplified. Consequently, to reduce that noise, a 5-point averaging method is used. 


2.1.2 Differentiation: 

After acquiring 5 samples, first differentiation YO and second differentiation Y1 are carried by using 
the equations (1) and (2). 

Yo [0] = X[0] — X[2] (1) 

Y, [0] = X[0] — 2xX[2] + x[4] (2) 

These differentiated outputs YO[O] and Y1[O] are undergoing for further smoothing by using the 


equations (3) and (4) which are stored in YsO and Ys1I respectively. Summing and scaling of first and second 
smoothed derivatives YsO and Ys1 are calculated by using the equation (5) and is stored in YT. 


Yoo = = (1¥o [0] + 21¥o[1]] + 1¥o[2]1) (3) 
Yor = = (IY LO] + 21%, (2]1 + 14 (21) (4) 
Yr = 1.3 Yoo + 1.1Y,, (5) 


2.1.3 Detecting R-—peaks in QRS complex: 

For detecting R—peak, we need to find the threshold value Yth which is calculated as the 80% of the 
maximum peak Yp. 

The physiological signals are not having same amplitude for all subjects. These vary from person to 
person in nature. In order to obtain the peak value, we have to adopt a new technique. [nitially, when the 
signal is interfaced to the microcontroller, the first five records are used to obtain the peak value. This peak 
value Yp is subsequently used for obtaining the threshold value. 

Yin = 0.81 * Y, (6) 


For any point, YT [i] that exceeds or equals threshold value Yth is compared. 
Yrlil = Yin (7) 
Subsequently, if eight succeeding values are greater than or equal to the threshold Yth, then the 
segment might be thought of containing part of QRS complex. Correspondingly, a pulse is generated at the 


digital output pin of the microcontroller. This method is simple and straightforward to implement when 
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compared with the other algorithms and produces a pulse that is proportional to every R-peak of QRS 
complex. 


3. RESULTS AND ANALYSIS 
We have tested our device on different subjects, using Storage oscilloscope. Here one channel is 


connected to ECG amplifier and another channel is connected to R-peak detector. Therefore, one can witness 
ECG and detected R—peaks simultaneously. By using this algorithm, we have detected pulse for every R— 
peak of QRS complex [12-13]. 

Below, we are illustrating ECG signals along with their detected peaks which are generated 
accurately. In Figure 4, the R—-R distance of subject 1 is 0.7 sec. The distance between the corresponding 


pulses is also 0.7 sec. 
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Figure 4. ECG signal in channel | of oscilloscope Figure 5. ECG signal in channel | of oscilloscope 
and pulses of corresponding R—peaks in channel 2 of — and pulses of corresponding R—peaks in channel 2 of 
subject 1. subject 2. 


From Figure 6, the R-R distance of subject 3 ECG is 0.75 sec and the measured distance between 
the corresponding pulses is also same which is 0.75 sec. 
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Figure 6. ECG signal in channel | of Oscilloscope and pulses of corresponding R—peaks in channel 2 of 
subject 3. 


Since the derivative of the signal is calculated as the rate of change of signal, sometimes it may give 
rise to two peaks which are close to each other. The second pulse can be rejected in the algorithm. 
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4. CONCLUSION 

In this paper, we have processed ECG signal from different subjects. In each case, detected pulses 
for each R—peak was found to be accurate. This proposed method can tolerate noisy environments and is 
simple to implement when compared with other methods such as template matching techniques, wavelet 
based techniques, digital filter based techniques and nonlinear transformation based techniques. The 
proposed device is portable, dynamic and can be implemented on low-end microcontroller. Therefore, we 
can find the R—-peaks instantly and economically. This can be used in many other applications such as heart 
rate, heart rate variability and in evaluating Pulse wave transit time. 
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